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; 1 INTRODUCTION 



We present new results on the gravitational lensing shear and magnification power 
spectra obtained from numerical simulations of a flat cosmology with a cosmological 
constant. These results are of considerable interest since both the shear and the magni- 
fication are observables. We find that the power spectrum in the convergence behaves 
as expected, but the magnification develops a shot-noise spectrum due to the effects 
of discrete, massive clusters and symptomatic of moderate lensing beyond the weak- 
lensing regime. We find that this behaviour can be suppressed by "clipping" of the 
largest projected clusters. Our results are compared with predictions from a Halo 
Model-inspired functional fit for the non-linear evolution of the matter field and show 
excellent agreement. We also study the higher-order moments of the convergence field 
and find a new scaling relationship with rcdshift. In particular, the statistic S3 is 
found to vary as z - 2 00 ±o o 8 (where z s is the source redshift) for the cosmology stud- 
ied, which makes corrections for different median redshifts in different observational 
surveys particularly simple to apply. 



Key words: Galaxies: clustering — Cosmology: gravitational lensing 
numerical — Large-scale structure of Universe 



Methods: 
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Knowing the distribution and evolution of the large-scale 
structure in the universe, together with the cosmological pa- 
rameters which describe it, are fundamental to obtaining a 
detailed understanding of the cosmology in which we live. 
Studies of the effects of weak gravitational lensing in the 
images of distant galaxies are extremely useful in providing 
this information. In particular, since the gravitational de- 
flections of light arise from variations in the gravitational 
potential along the light path, the deflections result from 
the underlying distribution of mass, usually considered to 
be in the form of dark matter. The lensing signal therefore 
contains information about the clustering of mass along the 
line-of-sight, rather than the clustering inferred from galaxy 
surveys which trace the luminous matter. 

Most obviously, weak lensing induces a correlated dis- 
tortion of galaxy images. The magnitude of the correlations 
depends on the density parameter, Q m , and the value of the 
vacuum energy density parameter, Qy, for the universe, as 
these parameters reflect both the amount of mass and the 
rate of evolution of structure. Consequently, the correlations 
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depend strongly on the redshifts of the lensed sources, as de- 
scribed by Jain & Seljak (1997) and Barber (2002). Recently 
a number of observational results have been reported for the 
so-called cosmic shear signal, which measures the variances 
in the shear on different angular scales. Bacon, Refregier 
& Ellis (2000), Kaiser, Wilson & Luppino (2000), Maoli 
et al. (2001), Van Waerbeke et al. (2000a, b), Wittman et 
al. (2000), Mellier et al. (2001), Rhodes, Refregier & Groth 
(2001), Van Waerbeke et al. (2001), Brown et al. (2002b), 
Bacon et al. (2002), Hoekstra, Yee & Gladders (2002), Hoek- 
stra, Yee, Gladders, Barrientos, Hall & Infante (2002) and 
Jarvis et al. (2002) have all measured the cosmic shear and 
found good agreement with theoretical predictions. 

In addition to shearing, weak gravitational lensing may 
cause a source at high redshift to become magnified or de- 
magnified as a result of the amount and distribution of mat- 
ter contained within the beam. The degree of magnification 
is strongly related to the convergence, which represents a 
projection of the density contrast, S(x), at position x along 
the line of sight, and which is proportional to Q m . Conse- 
quently, measurements of statistics for the convergence are 
able to provide cosmological constraints, and comparisons 
of the power spectrum of the convergence with theoretical 
predictions and numerical values serve to validate our the- 
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oretical and numerical models (see, e.g., Moessner & Jain, 
1998, and Jain, 2002). 

Of particular importance for interpreting weak lensing 
statistics is the fact that the scales of interest lie largely in 
the non-linear regime (see, e.g., Jain, Seljak & White, 2000). 
On these scales, the non-linear gravitational evolution in- 
troduces non-Gaussianity to the convergence distribution, 
and this signature becomes apparent in higher-order mo- 
ments, such as the skewness. In addition, the magnitude of 
the skewness values is very sensitive to the cosmology, so 
that measurements of higher-order statistics in the conver- 
gence may be used as discriminators of cosmology. 

In this work, we have obtained weak lensing statistics 
from cosmological ./V-body simulations using an algorithm 
described by Couchman, Barber & Thomas (1999) which 
computes the three-dimensional shear in the simulations. 
The code has been applied to cosmological simulations with 
fim = 0.3 and Q.y = 0.7; cosmologies of this type will be 
referred to as LCDM cosmologies. To obtain the required 
statistics on different angular scales, the computed shear 
values have been combined (using the appropriate angu- 
lar diameter distance factors and accounting for multiple 
deflections) along lines of sight arranged radially from the 
observer's position at redshift 2 = 0. Detailed results are 
presented for background sources at 14 different redshifts 
(z s = 0.1 to 3.6) and angular scales from 1' to 32'. 

As a test of the accuracy of non-linear fits to the con- 
vergence power we compare the numerically generated con- 
vergence power spectra with our own theoretically predicted 
convergence spectra based on a Halo Model fit to numerical 
simulations (Smith et al., 2002). We also investigate the sta- 
tistical properties of the magnification power spectrum and 
test predictions of the weak lensing regime. We also report 
on the expected redshift and scale dependence for higher- 
order statistics in the convergence. 

A brief outline of this paper is as follows. In Section 2, 
we define the shear, reduced shear, convergence and mag- 
nification in weak gravitational lensing and outline how the 
magnification and convergence values are obtained in prac- 
tice from observational data. In Section 3 we describe the 
relationships between the power spectra for the convergence, 
shear and magnification fluctuations, and how the power 
spectrum for the convergence relates to the matter power 
spectrum. We also describe our methods for computing the 
convergence power in the non-linear regime. Also in this 
Section, the higher-order moments of the non-linear con- 
vergence field are defined. The numerical procedure we use 
to generate the shear, convergence and magnification fields 
from the simulations are presented in Section 4, while in Sec- 
tion 5 we present our results for the numerical and theoreti- 
cal comparison of the convergence power spectra, the power 
in the magnification fluctuations, and the higher-order mo- 
ments, particularly the S3 statistic. Finally we discuss the 
results and present our conclusions in Section 6. 



2 WEAK LENSING FIELDS 
2.1 Weak shear 

Ellipticity measurements of observed galaxy images can be 
used to estimate the lensing shear signal. One definition for 



the ellipticity (see, e.g., Blandford et al., 1991, and Bartel- 
mann & Schneider, 2001) is the complex ellipticity, 



Q11 - Q22 + 2iQ 12 



Q11 + Q22 + 2(QnQ 22 - Q? 2 )2 



(1) 



in which Qij is the tensor of second brightness moments for 
a fixed isophotal contour. 

The "reduced shear," g, for a galaxy image at angular 
position 6, is defined by 



9(0) 



1(0) 
1 - k(0) 



(2) 



where 7 is the complex shear and k is the lensing con- 
vergence. Both 7 and n are obtained as projections from 
the values of the second derivatives of the lensing potential 
along the light path. Their detailed definitions and their re- 
lationships to the lensing potential are given by Schneider, 
Ehlers & Falco (1992) and summarised by Barber (2002). 
The transformation between the source ellipticities, e^' , and 
the image ellipticities, e, is given by 



(3) 



for I 9 l< 1- The asterisk in equation (3) denotes the complex 
conjugate. In the case of weak lensing, for which k, | 7 | and 
I g I are much less than unity, the transformation reduces to 



e M +9 



(4) 



for low intrinsic source ellipticites. 

The intrinsic ellipticities of given galaxies are not 
known, but averaging over the binned galaxy distribution, 
and assuming random ellipticities, yields a net lens shear: 



7 : 



<£>■ 



(5) 



This equality suggests that for weak lensing the variances 
in both the shear and the reduced shear for a given angular 
scale are expected to be similar. However, from numerical 
simulations, Barber (2002) has given explicit expressions for 
both as functions of redshift and angular scale, which show 
the expected differences. 

It is also possible to reconstruct the convergence from 
the shape information alone, up to an arbitrary constant, 
using methods such as those described by Kaiser & Squires 
(1993) and Seitz & Schneider (1996) for the two-dimensional 
reconstruction of cluster masses. Kaiser (1995) generalised 
the method for applications beyond the linear regime. 

Drawbacks to the reconstruction method arise from con- 
tamination by intrinsic galaxy alignments (Pen, Lee & Sel- 
jak, 2000, Brown et al., 2002a, Crittenden et al., 2001, Cate- 
lan, Kamionkowski & Blandford, 2001, Mackey, White & 
Kamionkowski, 2002, Heavens, Refregier & Heymans, 2000, 
and Croft & Metzler, 2001), although these can be statisti- 
cally removed if redshift information is available ( Heyman 
& Heavens, 2002, and King & Schneider, 2002). In addition, 
map-making of the convergence field over a finite area suf- 
fers from non-local effects due to missing information beyond 
the survey area (Bacon & Taylor, 2002). For these reasons 
it is useful to have an alternative method for estimating the 
convergence. 
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2.2 The magnification effect 

The lensing magnification, /i, can be computed directly from 
1 



H = (Idct^ir 1 = 



|(l-K)2_ 7 a|' 



(6) 



where A is the two-dimensional Jacobian matrix which de- 
scribes the mapping of a source onto its image. The effect of 
magnification on galaxy source counts results in a decrease 
due to the increase in the lensed image area, and an increase 
in counts due to the brightening of galaxies allowing their 
inclusion in a flux-limited catalogue. For a power-law flux 
distribution, the effect of lensing is (Broadhurst, Taylor & 
Peacock, 1995, and Bartelmann & Schneider, 2001): 



n (> S,z)~ f i(z) a - 1 n(>S,z), 



(7) 



where n(> S,z) and no(> S,z) are the lensed and un- 
lensed number of galaxy images per unit solid angle with 
flux greater than S and with redshift within dz of z, and a 
is the power-law exponent of S. Hence with calibration of the 
underlying number count amplitude and slope it is possible 
to estimate the magnification averaged over redshift. 

Alternatively, magnification values may be obtained 
from the change in image sizes at fixed surface bright- 
ness. This method is described in detail by Bartelmann and 
Narayan (1995) and summarised concisely by Bartelmann 
and Schneider (2001). See also Jain (2002) for a recent dis- 
cussion. 

Estimates of the lensing magnification based on number 
counts suffer from noise arising from the intrinsic clustering 
of the source galaxies, if redshift information is not available 
(e.g., Broadhurst, Taylor & Peacock, 1995, and Bartelmann 
& Schneider, 2001). With redshift information one should 
either select galaxies at different redshifts to remove the in- 
trinsic clustering signal, or use the brightening effect behind 
structure (Dye et al., 2001). Size distortions may also be af- 
fected if size is correlated with the environment. Again this 
may be removed by selecting galaxies at different redshifts. 

While the magnification signal-to-noise ratio is gener- 
ally poorer than that of the shear reconstruction method, it 
is valuable as an independent signal. Furthermore, if good 
values for the convergence are available, then the higher- 
order statistics are potentially very fruitful in discriminating 
amongst cosmologies (see, e.g., Bernardeau, Van Waerbeke 
& Mellier, 1997, and Jain et al., 2000). 

Finally, the convergence can be expressed (see Jain et 
al., 2000, for example) as a projection of the density contrast, 
6(0,x 3 ): 



2k(0,x s 



dx3 D( X3 )D( Xs - X3 )S(e^) As) 



D(x s 



where is the direction angle of the source at distance x s 
along the coordinate direction X3, Ho is the Hubble param- 
eter, D(xs), D(x s — X3) and D(x s ) are the angular diameter 
distances to position X3, x s to x$, and to position x s , re- 
spectively, and a is the scale factor. 

Although the lens convergence depends on the distance 
to the source galaxy this dependence is usually lost by aver- 
aging over the source distribution. However with redshift in- 
formation the full three-dimensional distribution of the den- 
sity field can be fully recovered (Taylor, 2001, and Bacon & 
Taylor, 2002). 



3 STATISTICAL PROPERTIES OF LENSING 
FIELDS 

3.1 The power spectra 

The convergence field can be expanded in two-dimensional 
Fourier modes on a flat-sky: 



k(£) = J 



d 2 GK(e)e l££ , 



(9) 



where £ is the angular wavenumber. The two-point correla- 
tion of these modes defines the power spectrum, C\ 



{k(£)k{£')) = (2n) 2 CrS D (£ - £')■ 



(10) 



Since the convergence can be expressed in the form of a 
projection of the density contrast, the power spectrum for 
the effective convergence can be obtained in terms of the 
matter power spectrum, Ps(k), for the density contrast (e.g., 
Bartelmann & Schneider, 2001): 

f WH .... fW(w) N 1 



a 



4c 4 



dw 



a(w) 



Ps 



(11) 



where the integral is evaluated over the comoving radial co- 
ordinate, to, from to the horizon, w H , defined by 



w H 



f dz 



(12) 



where H(z) = Ho[(l + z)' i Q, m + S7v] 1,/2 and where we have 
assumed a spatially flat universe; the weighting function, W, 
can be expressed as 

™ „.,/ 



W(w) 



dw' G(w')- 



(13) 



where G(w)dw = p z (z)dz, where p z (z) is the matter distri- 
bution function. 

On angular scales smaller than about 10', the total 
power in the effective convergence per logarithmic interval 
in wavenumber, 1(1 + l)C" re /(27r), is dominated by galaxy 
clusters. Jain et al. (2000) show that for the scales of interest 
in weak lensing, £(£ + 1)C( K / '(2ir) lies almost entirely in this 
regime, and there is significant enhancement (approximately 
an order of magnitude) of the power over linear predictions 
on scales below I ~ 10 4 . 

The two-point statistical properties of the shear and 
convergence are closely related. In the flat sky approxima- 
tion the components 71 and 72 of the shear are related to 
the effective convergence, k, in Fourier space (e.g., Barber, 
2002): 



7~i 2 W+72 2 W =k 2 {£). 



(14) 



Then it is clear that the power spectra for the shear, CJ 1 , 
and the convergence, C( K , are the same in the case of weak 
lensing. 

Finally, in the weak lensing regime, equation (6) reduces 

to 



H = 1 + Sfj, ~ 1 + 2k, 
so that 

Cf* 1 ~4CT = 4C7 7 , 



(15) 



(16) 



where C* 1 * 1 is the power in the magnification fluctuation, 
Sfi. This relation is questionable, since in general it is not 
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obvious that the weak lensing regime will hold over the whole 
magnification field. In particular, the denominator of the 
expression for n in equation (6) is sensitive to the value of 
k near peaks in the convergence field. We shall test this 
dependence in Section 5.2. 



3.2 Non-linear evolution of the matter field 

So long as the gravitational potential field is small the per- 
turbations in the matter density field can be fully non-linear. 
The three-dimensional gravitational potential, 3>(cc), is re- 
lated to arbitrary perturbations in the mass-density field by 
Poisson's equation, 



s 3 (e) 



WO))* 



(18) 



(V 2 + 3U 2 Ok)$(x) = 



(17) 



where Qk is the curvature energy-density parameter. 

A number of useful techniques have emerged to allow 
us to transform from the linear matter power spectrum to a 
fully non-linear spectrum. The first of these originated from 
the work of Hamilton et al. (1991) who considered the evo- 
lution of the matter correlation function, and was extended 
by Peacock & Dodds (1996) to account for the non-linear 
evolution of the matter power spectrum. These methods are 
based on the conservation of mass and a rescaling of physi- 
cal lengths due to gravitational collapse. These fits generally 
are accurate to around the 10% level. 

More recently Peacock & Smith (2000) and Seljak 
(2000) have developed a model for the non-linear evolution 
of the power spectrum based on the random distribution 
of dark matter haloes, modulated by the large-scale mat- 
ter distribution. This Halo Model for non-linear evolution 
reproduces the matter power spectrum of N-body simula- 
tions over a wide range of scales and has the advantage of 
relating the linear and non-linear power at the same scale. 
Given the utility of the non-linear fitting formula, Smith 
et al. (2002) have presented a new set of fitting functions 
based on the Halo Model functional form and calibrated to 
a set of TV-body simulations. These prove to be far more 
accurate (~ 1%) than previous formula. However these fits 
are only as accurate as the underlying simulations used in 
the fitting, which in this case were provided by the VIRGO 
Consortium^. In this paper we shall use the latter, more ac- 
curate fits based on the Halo Model functional form. While 
Smith et al. (2002) have compared their fit to the non-linear 
matter power spectrum, here we compare for the first time 
the predicted convergence power with the results of our sim- 
ulations in Section 5.1. 



3.3 Higher-order moments in the convergence 

In addition to two-point statistics, higher-order statisti- 
cal properties of the lensing fields are also of interest as 
non-linear evolution of the density field will introduce non- 
Gaussianity (e.g., Jain et al., 2000). Bernardeau et al. (1997) 
have investigated analytically the dependence of higher- 
order moments in the convergence on the cosmological pa- 
rameters, and in particular have discussed the ratio 



The significance of the S3 statistic is that it is expected to 
be independent of the normalisation of the power spectrum, 
and can also be shown to be rather insensitive to the angular 
scale. 

In the case of Sly = 0, Bernardeau et al. (1997) have 
shown that 



(19) 



for z s — 1. The £l m dependence is slightly weaker for sources 
at high redshift, and at low redshift S3 becomes approxi- 
mately inversely proportional to Q. m . The redshift depen- 
dence of S3 in an Einstein-de Sitter cosmology is approxi- 
mately z' 1 - 35 . 

Jain et al. (2000) have investigated the values for S3 
in different cosmologies for sources at z s = 1, including an 
LCDM cosmology using TV-body simulations based on recon- 
structing the convergence values from the shear, and includ- 
ing the effects of noisy data. Using various statistics based 
on the reconstructed convergence, they show that there are 
clear differences in the S3 values between the LCDM cos- 
mology and an open cosmology, and also claim that Q. m can 
be constrained to within an uncertainty of 0.1 - 0.2 in a 
deep survey of several square degrees. We shall study the S3 
statistic in more detail in Section 5.3. 



4 NUMERICAL PROCEDURE 

To evaluate the weak lensing statistics, we have applied the 
algorithm for computing the shear in three dimensions, as 
described by Couchman et al. (1999) to the cosmological TV- 
body simulations of the Hydra Consortium^ produced us- 
ing the 'Hydra' TV-body hydrodynamics code (Couchman, 
Thomas & Pearce, 1995). Simulations of the LCDM Dark 
Matter only cosmology were used with Q m = 0.3, ilv = 0.7, 
power spectrum shape parameter T = 0.25 and normalisa- 
tion, as, on scales of 8ft _1 Mpc of 1.22. The number of parti- 
cles, each of mass 1.29 x lO 11 /?. -1 solar masses, was 86 3 and 
the minimum value of the (variable) particle softening was 
chosen to be 0.0007(1 + 2) in box units. The simulation vol- 
umes had comoving side dimensions of 100/i _1 Mpc. To avoid 
obvious structure correlations between adjacent boxes, each 
was arbitrarily translated, rotated (by multiples of 90°) and 
reflected about each coordinate axis, and in addition, each 
complete run was performed 10 times. 

The general procedure for establishing the locations 
within the simulations for the computations of the shear 
and for computing the values of the elements of the shear 
matrices, is as described by Barber (2002), with the multi- 
ple lens plane theory being applied along the lines of sight. 
In this work, a total of 455 x 455 lines of sight were used 
and 300 evaluation locations for the three-dimensional shear 
along each line of sight in each simulation volume, thereby 
allowing regular sampling of the 2.6° x 2.6° field of view. 
With this number of lines of sight, the angular resolution 
equates to the minimum value of the particle softening at 
the optimum redshift, z = 0.36, for lensing of sources at 



t http://star-www.dur.ac.uk/~frazerp/virgo/virgo.html 



t (http://hydra.mcmaster.ca/hydra/index.html) 

© 0000 RAS, MNRAS 000, 000-000 



Simulated power spectra 5 



a redshift of 1. To allow for the larger angular size of the 
minimum softening at low redshifts and also for the range 
of particle softening scales above the minimum value, a res- 
olution limit of 1 arcminute has been adopted for the data 
analyses. 

A total of 14 source redshift slices were selected to give 
good statistical coverage of the redshifts of interest. These 
were redshifts of z s = 0.10, 0.21, 0.29, 0.41, 0.49, 0.58, 0.72, 
0.82, 0.88, 0.99, 1.53, 1.97, 3.07 and 3.57, corresponding to 
the redshifts of the simulation boxes. Hereafter, we shall only 
quote these redshifts to one decimal place for brevity. 

From the statistics computed in each of the 10 simula- 
tion runs, we computed the variances, the skewnesses, the 
statistic S3, and the power spectra for each of the source red- 
shifts. The two-point and higher-order moments were com- 
puted on angular scales of l'.O, 2'.0, 4'.0, 8'.0, 16' .0 and 
32'. using a top-hat filter, and the power spectra values 
were computed for a set of 15 wavenumber bins, spaced log- 
arithmically. Here we report on the results for these statistics 
as computed for the effective convergence and the magnifi- 
cation fluctuation (i.e., the departure of the magnification 
value from unity). We computed these statistics from the 
convergence values directly, rather than on convergence val- 
ues reconstructed from the shear. 

For the power spectra, the square of the absolute values 
of the Fourier transform of the convergence or magnification 
fluctuation were normalised by multiplication by j^ji -2-k£ 2 , 
where L is the angular length of the map side in radians. 
With this definition of the power spectrum, the results are 
essentially equivalent to £(l + 1)C^ K /(2tt). 

The computed values for the required statistics from 
each of the N = 10 runs were averaged, and the errors on 
the means of la/\fN determined. However, since we found 
the power in the magnification fluctuation to be a very noisy 
statistic, we have determined the median values and com- 
puted the errors on the medians for 
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Figure 1. Simulated 1(1+ l)C^ re /(27r) values (full lines) together 
with the linear (dot-dashed lines) and non-linear (dashed lines) 
predictions based on the Smith et al. (2002) Halo Model-inspired 
fitting formula, for z s = 0.5 (lowest set of curves), 1.0 (middle 
set) and 2.0 (uppermost set). 

the variance in the convergence field to be greatest. In terms 
of the simulation data for redshifts greater than 1, our field 
of view is effectively smaller because we make use of the peri- 
odicity in the simulations when lines of sight extend outside 
the volumes. Consequently, we sample in real terms a more 
limited area of sky, which suggests that the quoted error 
bars for the high redshift data may be understated. For the 
Halo-model, we simply note that its predictions may not 
adequately account for the presence of voids and detailed 
structure in the dark matter halos which become increas- 
ingly important in terms of their gravitational lensing effects 
for high redshift sources. 



5 RESULTS 

5.1 Comparison of the convergence power spectra 
with theoretical predictions. 

We have verified from our simulation results the equiva- 
lence of the shear and convergence power spectra. The shear 
statistics, in particular the shear variance, and their redshift 
dependence are described fully in Barber (2002). Here we 
present our results for the convergence power spectra for 
sources at the different redshifts obtained from our weak 
lensing simulations and from the linear and non-linear the- 
ory predictions. 

The convergence power spectra values from the weak 
lensing simulations are available for all the selected source 
redshifts. Figure 1 shows a comparison of the power spectra 
from our simulations with linear theory and the non-linear 
convergence power based on the Halo Model-inspired fitting 
functions of Smith et al. (2002) for sources at redshifts 0.5, 
1.0 and 2.0. We see that our numerical results lie primarily 
in the non-linear regime and that there is good agreement 
between the non-linear predictions and our simulations. No- 
table discrepancies are apparent only for the highest red- 
shifts and in the highly non-linear regime, where we expect 



5.2 The power in the magnification fluctuations 

We saw in Section 3.1 that the power in the magnification 
fluctuation, Cf 15 ^, is expected to be four times the shear 
and the convergence power. This result follows from mak- 
ing the assumption of weak lensing, where 5fj, (= // — 1) is 
taken to be equal to 2k. At low redshift, where the lens- 
ing may be considered to be weak on most scales, Cf" 5M /4 
and C( K are consistent with each other and with the weak 
lensing approximation. However, as the source redshift in- 
creases, we see increasing departures from this equality. 
Figures 2, 3, 4, 5 and 6 compare £(£ + l)C^/(8n) and 
1(1 + l)Cf7(27r) for source redshifts of 0.5, 0.8, 1, 1.5 and 
2 respectively. The departures are seen at first just on small 
scales (where the lensing is expected to be stronger), until, 
finally, at the highest redshifts, the departure is not only sig- 
nificant, but becomes consistent with a Poisson distribution 
for the magnification fluctuations on all scales. In support of 
this effect, we noticed the presence of large variations in the 
magnification power from one simulation run to the next for 
sources at high redshift. Because of these large variations, 
the curves have been plotted for the median power values in 
the magnification, rather than the mean, which would have 
been strongly influenced by the large power values. 
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for z s = 0.8. 



/(8tt) and mean £(t+l)C£ K /(2tv) Figure 5. Median l{l + l)Cf"*7(87r) and mean i{i + l)C£ K /(2n) 



for z s 



1.5. 



The departures in Cf *74 from CT arise from medium 
and strong lensing events. By making cuts in the data at 
the extreme ends, we are able to show the effects of high 
magnification events. Figure 7 reproduces the magnifica- 
tion and convergence data from Figure 4 and also shows 
£(£ + l)Cf*/(87r) and £(£ + 1)C? K / (2ir) after making cuts 
in the data beyond 6<r and la in their distributions. When 
the extreme magnification (and convergence) events are ex- 
cluded, we see that the apparent Poisson distribution in the 
magnification is removed and the curve reverts to the more 
"normal" behaviour expected from lower redshift sources. 
Thus we can conclude that it is the very high magnifica- 
tions which dominate the signal, although there is clear ev- 
idence of departures even for moderate levels of magnifica- 
tion. Our simulations have a relatively high value for the 
normalisation, as, when compared with very recent deter- 
minations (e.g., Brown et al., 2002b) and we might expect 
this to lead to a higher number of clusters than in cosmolo- 
gies with a lower as- Consequently, the shot-noise effect we 



describe may be moderated, although still present, in a uni- 
verse with a lower normalisation value. 



5.3 Higher-order moments 

In Figure 8 we plot the Ss{9, z s ) statistic for the convergence, 
which is a known discriminator of cosmology and indepen- 
dent of the matter power spectrum normalisation, for four 
different source redshifts and for angular scales from l'.O to 
32'.0. 

The redshift variation of Ss(6,z a ) for the different an- 
gular scales is displayed in Figure 9. On scales larger than 
8'.0 the errors in the measurements are greater, and these 
have not been plotted for clarity. By writing 

S 3 (e,z s ) = a(6)z He \ 



with 6 expressed in arcminutes, we find 
a(8) = (116 ±6)- (7.4 ±1.9)6* 



(20) 



(21) 
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Figure 7. The three solid curves represent the median l(t + 
l)C^'" 5 ' l /(87r) for the full magnification distribution (upper curve), 
and the cut distributions beyond 6cr (middle solid-line curve) and 
1(7 (lower solid-line curve) for z s = 1. These are compared with 
the equivalent curves for the mean 1(1+ l)C" K /(27r) for the con- 
vergence (dottcd-linc curves). 



and 

b(6) = -(2.00 ±0.08), 



(22) 



i.e., b is independent of 6 at a constant value, so that we 
may write, without the error values for clarity, 

-2.00 



S 3 (0, 



(116- 7.40)2" 



(23) 



for 1' < 9 < 8'. The fitting to this formula is excellent for 
all source redshifts up to and including 1, appropriate to 



most recent galaxy surveys. The precise 



dependence 



for S3 is also particularly simple to apply when combining 
or making comparisons with different surveys in which the 
median redshifts for the galaxies may be different, making 
this relationship particularly valuable. 

In Figure 10 we have determined the values of 



Figure 9. S:s(9,z a ) vs. z s at angular scales of l'.O, 2'.0, 4'.0 and 
8'0. 



{S30 2 K )(8,z s ), where a 2 represents the variance in the con- 
vergence, for the different source redshifts. For a given angu- 
lar scale, this product was found to be almost independent 
of redshift for an Einstein-de Sitter universe by Bernardeau 
et al. (1997) and we find a similar behaviour in the LCDM 
cosmology here, particularly at high redshift. Because of the 
independence from redshift, there is no need to adjust this 
statistic when making comparisons amongst different sur- 
veys, provided the median redshifts are not too small. This 
fact makes this statistic particularly useful for the discrimi- 
nation of cosmologies. 

Our results for Ssa 2 at low redshift are supported by 
combining the shear variance, a 2 = {-y 2 (9,z a )) ~ a 2 , with 
the expression for £3(0, z a ) given by equation (23). Barber 
(2002) has shown that 

( 7 2 (0,, s )) = ( K 2 (0,, s ))oc 2 r 7±o - 04 , 



(24) 

for source redshifts z s < 1.6 and angular scales of 2'.0 < 6 < 
32'. 0, so that the combination 
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Source redshift, z s 



Figure 10. S^cr^ vs. z s for angular scales of l'.O, 2'.0, 4'.0 and 
8'.0. 

S 3 O*{0, Zs) OC z 2-07±0.04 2 -(2.00±0.08) (25) 

predicts the very slowly rising function of redshift for z s < 1 
and 2'.0 < 9 < 8'.0, which are the ranges of applicability 
common to both results. 



6 DISCUSSION AND CONCLUSIONS 
6.1 Discussion 

We have shown that the convergence power spectrum val- 
ues computed directly from the weak lensing simulations 
show remarkably good agreement with the non-linear pre- 
dictions for the convergence power based on the Smith et al. 
(2002) Halo Model-inspired fitting formulae. Whilst we have 
not compared the predictions for other cosmological mod- 
els in this paper, the same numerical data, together with 
equivalent data for an open cosmology, have been compared 
with the predictions from Hierarchical models for the density 
field. Valageas, Barber & Munshi (2003) have shown excel- 
lent agreement for the full probability distribution functions 
for the shear and Barber, Munshi & Valageas (2003) have 
shown agreement with predictions for the convergence and 
higher-order moments. Consequently, these several agree- 
ments give us considerable confidence in our lensing pro- 
cedure and the resulting weak lensing data, which we have 
here analysed in a number of ways. 

Our results for the power in the magnification fluctu- 
ation suggest that the observed magnifications should be 
treated with some care as an estimate of the convergence 
power. This is especially true on small scales and at high 
redshift, since the magnification is very sensitive to locally 
high values of the convergence (arising from the tail of the 
skewed convergence distribution) and may vary significantly 
across small angular intervals. It further confirms that the 
weak lensing regime for the magnification should be treated 
with care, even for galaxy surveys with a mean redshift of 
0.8. Worse, at a redshift of 1, where many surveys have been 
undertaken or planned, the magnification power may be ex- 
tremely noisy. Consequently, evaluation of the convergence 



from magnification measurements is clearly likely to be bi- 
ased by isolated high peaks in the convergence and shear 
fields. 

The origin of this effect is the appearance of a few 
regions of very high magnification in the simulations. 
The spectrum is dominated by the strongest magnification 
events, although the signal shows moderate departures from 
the expected values even for medium magnifications. Since 
in an LCDM universe clusters form earlier and have time 
to relax, we expect large clusters to generate a large shear 
field. It is these largest structures with the highest shears, 
at the tails of the shear distribution, which dominate the 
magnification. These large shear-producing structures are 
rare, and so can be assumed to have a Poisson distribu- 
tion. Hence the effect is to produce a shot-noise effect in the 
observed magnification power spectrum. Since the magnifi- 
cation is a non-linear function of the shear and convergence, 
this will not appear in the individual shear and convergence 
power spectra, but should be more apparent in their higher- 
order moments. Another consequence of this is that the cos- 
mic variance, from realisation to realisation, should fluctu- 
ate quite wildly, which is the case here, and motivated the 
use of median statistics, rather than means which are more 
susceptible to this effect. 

Since the magnification effect results from lensing from 
the largest structures and since the formation of structure 
evolves at different rates in different cosmologies, we would 
expect similar effects to be present in other cosmologies, but 
with differing degrees of severity. We have not studied the 
effect in other cosmologies but we might reasonably assume 
that the LCDM cosmology would display large effects be- 
cause of the early development of structure. 

We have shown that the effect can be damped by "clip- 
ping" the magnification field by removing the highest peaks. 
Removing the magnification values beyond G-a in the dis- 
tribution removes most of the effect, showing that indeed 
it is the highest peaks causing this effect. One may won- 
der if this is a purely numerical effect, as in the simula- 
tions there is a high sampling of lines of sight, including 
the high-magnification regions around clusters. In reality we 
may not expect to see such extreme effects as often since 
galaxy surveys only Poisson sample the magnification field. 
In addition, if the normalisation, as, were lower than our 
chosen value, the number of clusters would be expected to 
be smaller, again leading to fewer examples of extreme mag- 
nifications. However, even after clipping the magnification 
values beyond 6-a there is still a residual systematic. Only 
after clipping the values beyond 1-a in the distribution do 
we remove this effect significantly. 

This all suggests that estimating the convergence power 
spectrum from magnification should be handled with care, 
unless evaluated at low redshifts. In addition, the estima- 
tion of higher-order statistics from the magnification will 
be equally affected by the shot-noise from individual, mas- 
sive clusters which will also increase the effects of sampling 
variance on any large-scale measurement of the magnifica- 
tion field. Hence, while magnification remains a useful tool 
for probing the mass distributions in individual clusters, its 
value as a statistical probe over large scales may be compli- 
cated by these non-linear effects. 

The results of Barber et al. (2000) for S3 on large angu- 
lar scales for different cosmological models, and the consis- 
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tent results of Jain et al. (2000) clearly show the strong Q m 
dependence of S3. Consequently, provided the convergence 
field is obtained by reconstruction from the shear, the S3 
statistic may be used to good effect to discriminate different 
cosmologies. In determining S3 from their simulations, Jain 
et al. (2000) quote values for sources at redshift 1 only. It 
is reassuring to note that the magnitudes of their S3 values 
on the different angular scales, determined in quite a dif- 
ferent way from ours (by reconstruction of the convergence 
from the shear data), are in good agreement at the small- 
est scales where non-linear effects are expected to be most 
pronounced. As the angular scale increases our values fall 
slightly more rapidly with scale. 

Our expression for the angular scale and redshift depen- 
dence for S3 clearly demonstrates the necessity of knowing 
the redshift distribution of the sources. This is the same con- 
clusion reached by Barber (2002) who reported the angular 
scale and redshift dependences for the shear variance. 

However, now that the redshift dependence of S3 is es- 
tablished for our LCDM cosmolgy (S3 oc z~ 2 ' 00 ), values de- 
termined from different surveys can be adjusted easily to a 
fixed median source redshift, enabling the combination and 
comparison of data from the different surveys. The only pro- 
viso here is that the redshift dependence of S3 is likely to 
be cosmology dependent, since the magnitude of this statis- 
tic on a given angular scale and for a specific distribution 
of sources is known to be a discriminator of cosmologies. 
We have not yet investigated the redshift dependence of S3 
in other scenarios, although Barber et al. (2003) have re- 
ported on various aspects of S3 from simulations and semi- 
analytical predictions in two different cosmologies. We would 
also hope to investigate the higher-order statistics for the 
convergence smoothed with a compensated filter in a later 
work, as this will allow more direct comparison with ob- 
served shear data. 

The statistic S30- 2 . for a given angular scale is much 
less sensitive to the source redshift, as Bernardeau et al. 
(1997) also found for the Einstein-de Sitter cosmology. Fig- 
ure 10 shows it to be approximately independent of redshift 
for high redshift sources and, using the results for the red- 
shift dependence of the shear variance from Barber (2002) 
together with the redshift dependence of S3 reported here, 
the statistic is shown to be only a very slowly increasing 
function of z s for z s < 1. Consequently, S3C 2 may prove to 
be very useful for surveys in which the redshift distribution 
of the sources is uncertain. 

6.2 Conclusions 

For the convergence power spectrum, £(£+l)C^ re /(27r), com- 
puted numerically from the lensing statistics in our simula- 
tions, we find excellent agreement with the values computed 
using a halo model. We can therefore have considerable con- 
fidence in our weak lensing statistics when applied to the 
shear, convergence, magnification and higher-order statis- 
tics. In addition, with this consistency, theoretical determi- 
nations of lensing statistics on a wider range of scales than 
can be achieved with numerical simulations may be used 
with increasing confidence. 

Our results for one quarter the power in the magnifica- 
tion fluctuations, £(£ + l)Cf'* 1 /(87r), show that the magnifi- 
cation is susceptible to the effects of discrete massive clusters 



and large variations across small angular intervals. These ef- 
fects occur specifically beyond the weak lensing regime, and 
become apparent in the magnification even at low redshift 
and on small angular scales. Consequently, determination 
of the convergence field from magnification data should be 
treated with special attention. 

Our simple mathematical description for S3(d,z s ), 
showing it to be closely proportional to z^ 2 for the LCDM 
cosmolgy, makes it particularly simple to compare and com- 
bine the results from different surveys in which the median 
redshifts of the galaxies may be different. 

Finally, we found directly, and through combination 
with the mathematical expressions for the shear variance 
and S3, that the combined statistic S30- 2 . is only a very slowly 
increasing function of redshift for low redshift sources and 
approximately independent of redshift at high redshift, as 
has been found in the Einstein-de Sitter cosmology. This 
statistic, therefore, may be useful in comparing and com- 
bining the data from different surveys with different galaxy 
redshift distributions. 
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